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Abstract 

A simple quantum model explains the Levy-unstable distributions for in- 
dividual stock returns observed by ref.[l]. The probability density func- 
tion of the returns is written as the squared modulus of an amplitude. 
For short time intervals this amplitude is proportional to a Cauchy- 
distribution and satisfies the Schrodinger equation with a non-hermitian 
Hamiltonian. The observed power law tails of the return fluctuations 
imply that the "decay rate", 7(5) asymptotically is proportional to |q|, 
for large \q\. The wave number q, the Fourier-conjugate variable to the 
return x, is interpreted as a quantitative measure of "market sentiment" . 
On a time scale of less than a few weeks, the distribution of returns in 
this quantum model is shape stable and scales. The model quantitatively 
reproduces the observed cumulative distribution for the short-term nor- 
malized returns over 7 orders of magnitude without adjustable parame- 
ters. The return fluctuations over large time periods ultimately become 
Gaussian if 7(9 ~ 0) oc q 2 . The ansatz 7(g) = br \J vn? + q 2 is found to 
describe the positive part of the observed historic probability of normal- 
ized returns for time periods between T = 5 min and T ~ 4 years over 
more than 4 orders of magnitude in terms of one adjustable parameter 
st ~ mbr oc T. The Sharpe ratio of a stock in this model has a fi- 
nite limit as the investment horizon T — ► 0. Implications for short-term 
investments are discussed. 

1 Introduction 

The quest for a quantitative statistical description of stock prices began more 
than a century ago with Bachelier's thesis 2 in which he described price move- 
ments by a random walk. Although Bachelier's random walk has since been 
modified, the effort to model stock price movements by a stochastic process 
continues. Computerized historic financial records have recently made a high- 
precision statistical analysis of short-term returns possible. For time scales that 
are less than a few weeks it is increasingly difficult to reconcile this wealth of 
new information with the conjecture that price fluctuations are described by 
a stochastic process with independent and identically distributed increments 
{iid). 
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To estimate the probability density function (pdf), { x t )i f° r * ne return 
Xj. (the change in the logarithm of the price S« of the stock (i) over a time 
period T), 

4? := ln[^ T /5 t W ] = ln[sg> r ] - ln[S«] , (1) 

the frequency of such changes was analyzed in ref.[l], using historical records 
that typically cover several years. Reducing the time horizon T to minutes 
offered the advantage of suppressing the effects of a changing macro-economic 
environment while at the same time improving the statistics. 



T=5min- Iday 
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Fig. 1. Fits to the observed average cumulative distribution for short-term normalized 
returns. The historic probabilities for time periods of T — 5 minutes (dots) and 
T = 1 day ~ 390 minutes (squares) are from the analysis of ref. [1] . Only the positive 
tails of the slightly skew distributions are shown. The lines correspond to cumulative 
distributions of the form given in Eq. (13) for three different ratios, r — (a/b) = 0, 1, 3, 
of the parameters a and b. 

The analysis of ref. [1] revealed that the distribution of returns for different 
stocks and for different time intervals T < 2 weeks are all shape-similar and 

(i) 

exhibit a somewhat unexpected power law fall off. More precisely: if ii T 

denotes the mean and Vj) the standard deviation of the returns over 
a period T for the stock (i), the cumulative distribution of the normalized 
returns g, 

N®(g) = E[4> - $ > vfg\ , (2) 
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was found 1 to not depend on the stock (i) nor on the time horizon T for 

(i) . . (i) 

5 min < T < 2 weeks. The mean, fi T , and standard deviation, Vj.\ thus were 
found to be the only statistically significant characteristics of the historical 
distribution of returns of any individual stock on time scales up to a few weeks. 
Scaling in the following denotes the fact that the distributions of the normalized 
returns are very similar. Although this would also be true for a random walk or 
any other Levy stable process, it is somewhat unexpected that non-Gaussian 
return distributions with a finite variance are shape-similar over vastly different 
time horizons T, for very different companies and in very different economic 
environments. 

The empirically observed scaling of the distributions permits a significant 
improvement of the statistics by averaging over the distributions of the nor- 
malized returns of many individual stocks 1 , 



One thus can reliably estimate the probability for return fluctuations that are 
about 100 times larger than average. 

The shortest time horizon investigated by ref.[l] is T = 5 min. The anal- 
ysis in this case follows the stocks of the 1000 US companies with the largest 
market capitalization over a 2-year period from January 1994 to December 
1995. The average cumulative probabilities for the normalized 5-minute re- 
turns found by ref.[l] is reproduced in the log- log plot of Fig. 1. Also shown is 
the average cumulative distribution of the normalized daily returns. The lat- 
ter were extracted 1 from the records of stock prices for about 16,000 individual 
companies (binned by market capitalization) over the entire 35-year interval 



The interested reader is referred to the original analysis in ref.[l] for fur- 
ther details. But I would like to emphasize that the observed scaling of the 
distributions of individual companies (and of groups of companies with differ- 
ent market capitalization) is crucial to the interpretation of the analysis. The 
average normalized cumulative distribution of Eq. (3) otherwise would differ 
qualitatively from the normalized distribution of any individual stock. With- 
out scaling, the error in the mean normalized distribution of Eq. (3) would be 
large and the data base of about 4 x 10 4 events per company would not allow 
any conclusions about the frequency of events that occur only a few times out 



Fig. 1 shows that the averaged distributions for different time intervals 
T < 1 day also scale extremely well. It is remarkable that normalized data sets 




# stocks 



(3) 
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for time intervals that differ by a factor of about 78 ~ (1 trading day)/5 min 
and involve different companies in different economic periods should show no 
statistically significant difference. 

Apart from this shape-similarity over vastly different time scales and many 
different companies, the most striking feature of the empirical probability dis- 
tribution is the pronounced power law over ~ 5 orders of magnitude, 

N T {g) oc g~ a , for 2 < g < fOO and T < f week . (4) 

The exponent a empirically is close to 3. For T = 5 min, the best estimate for 
a is 3.10 ± 0.03 for the positive tail and 2.84 ± 0.12 for the negative tail of the 
observed distribution 1 . For short time intervals T < 2 weeks the distributions 
are only slightly skewed, but this feature becomes more pronounced in the 
data with increasing T. The observed asymmetry of the distributions could be 
caused by a number of factors, such as the discreteness of prices or bankruptcy 
regulations, that all tend to mainly distort its negative tail. We do not model 
such effects and for the purpose of this investigation replace the empirical 
distribution by a symmetric one with the same positive part. 

The empirical estimates of the power law exponent are well outside the 
region for Levy-stable distributions, which would require 3 that < a < 2. 
The stability of the power law tails with an exponent a ~ 3 refute Bachelier's 
conjecture 2 and its generalizations. The observed distributions cannot result 
from an iid process. The pdf for T = 390 min ~ 1 day otherwise would have 
to differ considerably in shape from the pdf for T = 5 min; the 390/5 = 78- 
fold convolution of the latter is rather close to a Gaussian pdf and not at all 
shape-similar to itself. 

The observed auto-correlation time is of the order of a few minute^only. 
Memory effects that could explain the persistence of power law tails over days 
therefore probably are of higher order. Although GARCH-processes 4 in princi- 
ple model such correlations, it is quite difficult to obtain a shape-stable tempo- 
ral evolution with the observed power law in this manner 5 . A GARCH-process 
furthermore depends on parameters that determine the shape of the marginal 
distribution. In order to reproduce the observed scaling, one would have to 
contend that these shape degrees of freedom somehow are strongly correlated 
and very similar for different companies, macro-economic conditions, etc. 

Power law tails with exponents a ~ 3 have also been observed for the cu- 
mulative distributions of market indices 6 and of commodity prices 7 . Although 
the existence of fat tails was recognized early on, an economic explanation of 
the power law with an exponent a ~ 3 has only recently been proposed 3 . In 
this micro-economic analysis, the asymptotic power law reflects the distribu- 
tion and trading behavior of the largest investors in a stock, such as mutual- 
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and/or pension- funds. However, such an asymptotic analysis cannot, by it- 
self, explain the actual magnitude of the tails nor the observed scaling of the 
distributions. 

If one ignores the problem of temporal stability, the distribution of returns 
on a stock (including the power law tails) for any fixed (short) time interval 
T can be modelled by a subordinated stochastic process. Clark 9 originally 
proposed a subordinated stochastic process in his attempt to explain daily 
returns on cotton futures. He contended that such a process would give a 
better description of cotton futures than Mandelbrot and Taylor's suggestion 10 
of a Levy stable process. Only a limited amount of data (~ 2000 prices) 
on the daily closing of cotton futures was analyzed and the issue of shape- 
stability of the distribution of returns over longer (or shorter) time intervals 
than 1 day was not addressed. The daily returns on cotton futures indeed are 
better reproduced by a subordinated stochastic process than by a Levy stable 
one, mainly because the latter can give a leptokurtic distribution with finite 
variance. 

More recently, the Variance Gamma process 11 , a subordinated stochastic 
process with finite kurtosis, has been proposed to value options. A subordi- 
nated stochastic process with finite variance but infinite kurtosis that repro- 
duces the observed distribution of returns on a stock for any fixed time interval 
T < 2 weeks is constructed in the appendix. However, this iid process with 
finite variance converges rapidly to its Gaussian fixed point and does not ex- 
plain the shape-stability of the observed distribution on different time scales 
T. 

In short, the observed price fluctuations pose at least three theoretical 
challenges: 

1. To explain the scaling of the distribution of short-term returns of many 
individual stocks over ~ 7 orders of magnitude. 

2. To explain the power law behavior with a <~ 3 of the tails of the cumu- 
lative distribution for the returns. 

3. To explain the apparent temporal shape-stability of the distribution for 
the returns over time horizons from a few minutes to a few weeks and the 
slow convergence to a Gaussian distribution on much longer time scales. 

The functional form of the observed normalized distribution of short-term 
returns is accurately reproduced by the ansatz of the next section. In section 3 
the temporal evolution of the amplitude (the square root of the pdf) is found to 
be naturally shape-stable in a quantum model and the linear evolution operator 
(effective Hamiltonian) that reproduces the observed temporal evolution for 
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T < 1 day is obtained. Section 4 discusses and interprets a slightly modified 
effective Hamiltonian. It semi-quantitatively describes the temporal evolution 
of the distribution of returns over time scales up to four years. Section 5 
explores some of the implications of this quantum mechanical description and 
the results are summarized in section 6. 



2 Power Law Tails with a = 3 in Quantum Finance 

The recently proposed quantum description of financial markets 12 offers a sur- 
prisingly simple and transparent explanation for the observed distributions. 
The exponent a ~ 3 for the power law tails, in particular, is quite natural in 
this framework and is temporally stable. 

a = 3 in the cumulative distribution of Eq. (2) implies the asymptotic 
behavior, 

d 

p T {x 2 ~ oo) = - —N T (x/v T ) 



ccx, (5) 



of the pdf of a representative stock's return. The essence of a quantum de- 
scription is that the pdf is interpreted as the squared magnitude of a (possibly 
complex) amplitude <f>r(x), 

Pt(x) = |Mx)| 2 • (6) 

If the temporal evolution (and thus the pdf) is invariant under the trans- 
formation# x — \i — > \i — x, the amplitude either is symmetric or antisymmetric 
in x — /j. Since the probability of the mean return, pt(p-), does not vanish, 
4>t(x) in this case is symmetric and a function of (x — fi) 2 only. The asymptotic 
behavior of the pdf in Eq. (5) implies 

4>t{x 2 ~ oo) oc x~ 2 , (7) 

for the asymptotic behavior of the amplitude. Since <Pt{x) falls off like a power 
law for large values of x with an exponent that is (close to) a negative integer, 
we model the amplitude by a rational function, i.e. a function that is analytic in 
the whole complex plane apart from a finite number of poles. The position and 
strength of the poles in this case constitute the set of parameters that describe 
the pdf of an individual stock. The observed approximate scaling implies 
that the distribution of an individual stock's returns to first approximation is 
specified by its variance and mean. Our ansatz therefore should not involve 
too many poles. Every additional parameter is related to additional* statistical 

#This is an (approximate) " parity" -symmetry for reflection about the mean fi. 
*With more than two poles, the pdf can, for instance, be skewed. 
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properties that in general will break the scaling and distinguish between the 
shapes of the distribution of returns of individual stocks. An amplitude with 
a single pole corresponds to a pdf with infinite variance. The observed finite 
variance thus requires an amplitude with at least two poles and that the sum 
of residues must vanish. The amplitude then asymptotically falls off at least as 
fast as in Eq. (7). Note that only asymptotic power laws with an exponent a 
that is an odd integer can be modelled in this fashion and that an odd exponent 
a > 3 would require more than two polcs^. 

The simplest rational amplitude corresponding to a pdf of finite vari- 
ance therefore has two poles of opposite strength. It automatically has the 
asymptotic behavior of Eq. (7). The requirement that the amplitude is a sym- 
metric function of x — fj, further constrains the position of the two poles to 
z = fj, ± (a + ib). One thus is led to consider the ansatz, 

for the amplitude. The real parameters zt, a and b > generally will depend 
on the stock, economic era and the time interval T. Normalizing the pdf 
determines the constant Af up to an irrelevant phase, 

\M\ 2 = (a 2 + b 2 )- . (9) 

Note that the pdf is not normalizablc at b = 0. The normalized pdf corre- 
sponding to the ansatz of Eq. (8) is, 

2b(a 2 + b 2 ) 

PT{X > ~ 7r(((x- M ) 2 + 6 2 -a 2 ) 2 +4a 2 6 2 )' ( ' 

The three real parameters, /i, a and b > on which it depends are related to 
the mean, variance and the curvature at the mode of the density. The variance 
of Pt{x) is just the square of the distance of the complex poles from the mean, 

v 2 = \z-v\ 2 = a 2 + b 2 . (11) 

The second (independent) parameter is the ratio r — a/b, or, cquivalcntly, the 
phase of z — jj,. The curvature of the pdf at x = fx is 



92 I ^ 



(12) 



TFor a's that are not odd integers, the amplitude would have branch cuts 
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p is negative for a 2 < b 2 and the pdf in this case peaks at x — p. For a 2 > b 2 
the pdf is double-humped with two distinct maxima at x max = p ± \J a 2 — b 2 . 
pv 3 is a function of the ratio r = a/b only. 

The distribution for the normalized returns that corresponds to the pdf 
of Eq. (10) is, 



N T (g) = ± 



7r + arctan I — - — - I — arctan 




a+vg i 



(13) 

With Eq. (11), N T (g) is seen to depend on the shape parameter r 2 = (a/b) 2 
only. For large and small values of g, N T (g) has the expansions, 



Nt(q ~ oo) <~ ^—Q 3 



1 + ^ + 0(9 -., 



1 9h mi 3 

N T (g ~ 0) ~ - - — g - ^-g 3 + 0(g 5 ) . (14) 

Z 1TV 6 

[The relative correction to an asymptotic power law oc g~ 3 is never much more 
than 6/(5g 2 ) (or less than 1% for g > 11)]. The leading correction vanishes 
altogether for p — 0, or equivalently, r = 1. 

Fig. 1 shows distributions Nr{g) of Eq. (13) for ratios r — 0,1, 3, together 
with the empirical data. Apart from the overall change in the normalization of 
the tails, the parameter r qualitatively changes the shape of the distribution 
of normalized returns near g = 0. Note that the theoretical distributions, by 
construction, all have the same power law tails with an exponent a = 3, but 
that the normalization of the tails of the distribution for r = 3 is off by a 
factor of ~ 20. Any asymptotic analysis that predicts only the exponent?, by 
itself, therefore cannot provide an explanation of the observed distribution of 
returns. 

There is no convincing reason why parameters like r should be the same for 
different companies and/or economic eras or different time intervals T . Large 
variations in r would break the observed shape similarity of the distributions. 
Within the accuracy of the analysis of rcf. [1], the ratio r = ar/bx does not 
appear to depend strongly on T and should be similar for most stocks. The 
empirical data prefers small r < 1 and is consistent with r — 0, i.e. with an 
amplitude whose two poles have similar real parts. In this case the amplitude 
describing a mean return of pr with standard deviation vt = &t has two 
complex conjugate poles of opposite strength located at zt = Pt + ivr and 
z?p = pt — ivr- This is the minimal number of independent parameters the 
distribution of returns of an individual stock can depend upon. The normalized 
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probability distribution corresponding to this special case is the solid line in 
Fig 1, a student t-distribution with 3 degrees of freedom, 



1 



1 



2 



2v T 



Pt{x) = — 



2tt x — z t x — z T 



7^-/^)2+4)2 



(15) 



It fits the data exceedingly well. For T < 2 weeks, the returns of a stock 
in this sense are statistically described by the path the pole zt takes in the 
complex half-plane with Im zt > 0. The path begins at the origin zt=o = 
and long term investments are encouraged if the standard deviation vt is a 
concave function of the mean return /zt. 

3 Quantum Stability 

The good description of the empirical data by an ansatz for the pdf that is the 
square of the modulus of a rational amplitude could be considered fortuitous 
and in itself does not require a quantum model. But a quantum model of the 
dynamics does explain the apparent stability of the shape of this distribution 
over vastly different time intervals without the need to explicitly model memory 
effects. It also casts some doubt on the notion that stock market fluctuations 
are predictable if they do not follow an iid process. 

Quantum dynamics primarily describes the temporal evolution of the am- 
plitude - the evolution of the corresponding pdf follows from the relation 
in Eq. (6) and is non-linear. Stochastic processes that give Levy-unstable 
distributions of the type Eq. (10) (with the correct power law tails) over an 
extended period of time 5 ' 13 ' 14 tend to be arbitrary in the sense that other dis- 
tributions may almost equally well have been obtained. They either depend 
on quite a few carefully adjusted parameters 5,13 or are directed by a process 
whose fractal dimension would have to be a time independent characteristic 
of financial markets 14 . Unlike for physical phenomena, such "hidden" charac- 
teristics of financial markets are expected to slowly change over time and the 
distributions thus would depend on the economic era. The ansatz of Eq. (8) 
for the amplitude on the other hand turns out to be absolutely shape stable 
under a relatively simple and quite natural quantum dynamics. 

To better see this, consider the Fourier-transform of the amplitude of Eq. (8) 
with the normalization Eq. (9). Contour integration about the simple poles of 
the ansatz at z = fir ± («t + ibr) with b T > is elementary and gives the 
Fourier-transform 4>t{q), 




(16) 
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For very short times, when the mean and variance of the distribution are small, 
4>t(q) essentially is constant over a wide range of wave numbers q. In this case 
the corresponding amplitude Eq. (8) for the return is well localized about 
x = fi, i.e. the expected log-price. The variance of the pdf grows with the time 
horizon if ax and/or br grow in magnitude. As noted before, the cumulative 
distribution of the normalized returns remains strictly the same only if the 
parameters ar and br grow proportionally and the ratio r = (ot/&t) does not 
depend on the time interval T. 

The temporal evolution of a state \tp) in the Hilbert space of a quantum 
model is generated by some effective Hamiltonian operator H e a(t), 

\<p) T = Te- i fo dtii ^ t) \<p) , (17) 

where the symbol T denotes time ordering of the exponential factors (redun- 
dant if, as in the present case, the Hamiltonian operators for different times 
commute). To conserve probability in a complete Hilbert space, H e g(t) neces- 
sarily would have to be hermitian and states would evolve by a unitary rotation. 
However, if H e &(t) describes only the evolution in a subspacc of the Hilbert 
space - such as that spanned by a single share - probability can be "lost" to 
the complementary part of the full Hilbert space. H c fi{t) in this case is not 
hermitian. 

Comparing Eq. (17) with Eq. (16) suggests that the effective Hamiltonian 
for the time evolution of a share is diagonal in the Fourier-conjugate basis of 
| q) -states, with matrix elements, 

(q\H eS (t)\q') = 2tt% - q')(-qfi t - \q\a t - i\q\b t ) , (18) 

where the dot is shorthand for the derivative with respect to time, i.e. fit — 
d^t/dt, etc. 

In the language of particle physics, <j)r{q) oc {q\(p)r is proportional to the 
wave- function of a (massless) particle whose frequency ui{q) and decay rate 7(g) 
are both proportional to the magnitude of its wave-number q. If r = a/b does 
not depend on time, one can interpret &t as proportional to the " proper" time 
of a stock, i.e. as the monotonically increasing parameter that characterizes 
the evolution of the stock's price distribution, bj- need not be proportional to 
physical time or even to trading time. The proper time of a share could be 
proportional to some monotonically increasing quantity that is relevant to the 
actual trading dynamics of the stock, such as the overall number of trades or 
the number of traded shares. We will further examine this issue in section 5. 

From the available data one in fact cannot conclude that the frequency 
to(q) and decay rate j(q) of a share are proportional. As Fig. 1 shows, the 



10 



data implies only that the term proportional to is small compared to br\q\- 
Since ot and bx are the strengths of two commuting terms of the evolution 
operator, the contribution u>(q)ciT may simply be negligible, ar — is a trivial 
fixed point of the evolution that reproduces the empirical data rather well. 
[Note that this would not be the case if one could not separate the effects of 
the two parameters in H c g.] Assuming \ar\ <C &t in this sense is a "natural" 
approximation to the evolution that does not require any fine tuning. 

The time dependent proportionality constant ^/2nbT in the amplitude 
4>t(q) of Eq. (16) ensures that the pdf remains normalized at all times when 
&t > 0. It docs not appear in Eq. (17) and the norm of the state \<p)t de- 
cays with T if bx increases with time. The time dependent normalization 
factor arises because the effective "one-particle" Hamiltonian H e s describes 
the temporal evolution of a particular share of a stock. We are therefore com- 
puting the probability of a return x, if a certain share traded at time t is 
again traded at time t + T > t. The probability that any particular share is 
again traded (for any price) after a time interval T is just the normalization 
T ((p\ip} T = (2i:b T y 1 (ip\ip} of the state 8 \<p) T in Eq. (17). The probability 
that a particular share is traded after a time interval T thus decreases like 
1/&T, but the probability for a particular share's return under the condition 
that it is again traded after a time T of course remains normalized for all T. 
[The calculation is analogous to that of the probability of a decaying particle's 
change in position if it is observed after a time T] 

4 Ultimate Convergence to a Gaussian Fixed Point 

The exponent a ~ 3 and the apparent temporal shape stability of the distri- 
bution thus are linked by the fact that quantum theory primarily describes 
the temporal evolution of the amplitude, rather than of the transition prob- 
ability. Although the evolution of an amplitude of Cauchy type such as the 
one in Eq. (8) is shape stable, there are other possible fixed points. Of some 
interest is the Gaussian one, because one does expect the pdf to eventually ap- 
proach a Gaussian distribution. Note that the Fourier-transform of a Gaussian 
amplitude is a Gaussian and that a Gaussian amplitude implies a Gaussian pdf. 

The pdf pr{x) of individual stocks empirically resembles a Gaussian dis- 
tribution only after several years. The results of the analysis of ref.[l] for 
T = 16 days to T = 1024 days, or up to approximately 4 trading years, are 
shown in Fig. 2. The slow rate of convergence to a Gaussian fixed point sug- 

$ We will sec in section 5 that it is conceptually more appropriate to include a factor 2-7T& in 
the relation of Eq. (6) between the conditional probability density and the square of the 
amplitude, instead of normalizing the square of the amplitude. 
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gests that it may be due to macro-economic and other factors such as stock 
splits, buy-backs, etc., that are irrelevant for the short-term dynamics of a 
stock. Note that the positive and negative tails of the cumulative distribution 
do not approach the Gaussian fixed point at the same rate 1 . Although the 
approximation by a symmetric distribution becomes questionable on long time 
scales, it nevertheless should be possible to consistently describe the crossover 
from one (almost) fixed point to another within a single effective quantum 
model. For simplicity and because it is a fixed point and consistent with the 
short-term data, we consider only the case ax = 0, i.e. a stock whose evolution 
is completely dominated by its "decay rate" 7(g). 

T= 16- 1024 days 

I 0.1 

X 

1 0.01 

2 0.001 
I 0.0001 

1 

0.00001 

1 1.5 2 3 5 7 10 15 
Normalized return 




Fig. 2. The cumulative distribution for the normalized returns over time intervals 
from T — 16 to T = 1024 trading days. The empirical probabilities are from the 
analysis of ref.[l]. The full lines are distributions that correspond to the extended 
ansatz of Eq. (20) that best fit the data. They depend on the value of the parameter 
s = bra only. 

For long time horizons, the variance of the pdf and thus of the amplitude 
are large. The Fourier-conjugate amplitude 4>t(q) therefore is concentrated 
about q ^ 0, and the approach to a Gaussian fixed point can be ensured if the 
function 70(9) = \q\ in Eq. (16) is modified to one where the cusp at q = is 
replaced by a smooth quadratic dependence on q. It was argued previously 12 , 
that the effective decay rate 7(g) should be a quadratic function of q in the 
long wavelength limit, q ~ 0. The proposed modification for small values of q 
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will not affect the short-term power law behavior over a wide range of x if -f(q) 
for \q\ ~ oo approaches the function \q\ sufhciently rapidly. Replacing \q\ by 
the upper branch of a hyperbola with these asymptotes thus is a possibility 

M -> lm(q) oc \Jq 2 + to 2 - m . (19) 

The parameter to controls the transition between the Gaussian and power 
law regimes. For m = one recovers the previous case, whereas the density 
essentially is Gaussian for all but the shortest times when to is large. For m > 
the ultimate fixed point thus is Gaussian, but the convergence to it is slow for 
small to. The particular choice for the interpolating function in Eq. (19) is ad 
hoc, but a square-root dependence of the decay rate on the wave-number can be 
the result of diagonalizing a (anti-)hermitian 2x2 matrix - such square-roots 
are quite generic for quantum systems near level crossings. 




-7.5 -5 -2.5 2.5 5 7.5 
Number imbalance 

Fig. 3. The conditional expectation < N >$ of the number of trades for a given 
number imbalance <F The data points are reproduced from Fig. 3a) of the analysis 
of ref. [15]. The solid line represents the function \/& 2 + 0.49 — 0.7. 

Empirical support for the particular form of j m (q) in Eq. (19) is pro- 
vided by the analysis in ref. [15] of how the mean number of trades (iV)$ T in 
a time interval of T = 15 min depends on the imbalance between the buyer- 
and seller-initiated trades. Whether a trade is buyer- or seller- initiated was 
determined 15 from the prevailing quote shortly before execution of the trade 
using the procedure of C. Lee and M. Ready 16 ; the imbalance <&t then is the 
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difference in the number of buyer- and seller- initiated trades in the time in- 
terval. The data points in Fig. 3 are the results of ref. [15] from the analysis 
of the trading records of the 116 most frequently traded US stocks during the 
2-year period 1994-1995. 

If one assumes that the average volume of shares per trade does not depend 
strongly on <&t, the observations of ref. [15] can be interpreted and understood 
within the quantum model. Because the time interval of T = 15 min is rela- 
tively short, < N >$ T in this case is proportional to the "decay" rate 7 m (?) 
of a share in state \q). The observations in Fig. 3 indeed appear to lie on a 
symmetric hyperbolic curve. Furthermore, there appears to be little trading 
activity at vanishing imbalance <!>t ~ 0. It therefore is tempting to assume 
that the wave number q and the imbalance <&t are roughly proportional. This 
interpretation would associate the market dynamics captured by the imbalance 
$x with the Fourier conjugate variable to the return xt (up to a proportion- 
ality constant). It has potentially interesting and observable consequences in 
a quantum model. The pdf of 3>t in this case is related to the pdf of xt- the 
corresponding amplitudes are Fourier conjugate to each other (up to a scale). 
The amplitude </>t(<?) of Eq. (16) implies that for time intervals T < 2 weeks 
the normalized distribution of $t is E[$ > gv& T > 0] <~ ie~ 9 /^, where 
t>$ T oc 1/br is the standard deviation of $t- The model thus predicts that 
the probability of a large imbalance in the trades falls off exponentially, rather 
than like a power law. [Note that u$ T oc 1/bx is one version of Heisenberg's 
uncertainty relation. It here implies that periods with large fluctuations in the 
imbalance should show relatively small fluctuations in the returns, and vice 
versa.] 

Returning to the temporal evolution of the distribution of the returns, the 
modification of Eq. (19) leads to an amplitude of the form & , 



1 [°°dq , mb gi(myg+F) 
(j) T {x;m) = — / — cos(gx)exp[-&V<T + ™ ] = — T7 ; 9 , 9 , 

(20) 

where K\(z) is the modified Bcssel function of the third kind with the asymp- 
totic behavior 

K x {z~Q)~\/z, K x (z ~ oo) ~ 4/ ^-e~ z . (21) 



& We assume that aj- <C bj< and neglect any dependence on the parameter ay. This approx- 
imation is consistent with the short-term results. 
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Normalization of the transition probability requires that, 

|AA| 2 = r d A e -^^FT^ = m Kl{2mb) . (22) 
Jo n n 

For a mean return fi, the pdf for the return on a stock in this model then 
becomes, 



, h s U ( mb 2 K 2 (m^(x-^ 2 + b 2 ) 
p T (x;»,b,m) = \Mx-*rn)\ = vKl{2nh) {x _ „ )a + b 2 . ( 2 3) 

and depends on three, possibly T dependent, parameters /xt, br and tot- The 
variance of the returns is 

2 _ b2 f°° g 2rf g 26^/^+^ f24 x 

VT - mK^mb) J qi+m^' ' (24) 

For long time horizons (raft 3> 1), the limiting density is Gaussian with a 
variance v 2 (mb 1) ~ 6/ (2m). For short time intervals (to6 <C 1) the variance 
to first approximation does not depend on to, v 2 (mb <C 1) = b 2 (1 — 7tto6 + . . . . 

Since the variance of normalized returns is unity by definition, the cu- 
mulative distributions for the normalized returns depends on the parameter 
combination s = mb only. As shown in Fig. 2, this extended ansatz qualita- 
tively and to some extent even quantitatively reproduces the observed 1 positive 
tails of the (average) cumulative distributions of the normalized returns. The 
distributions for time periods of T = 16, 64, 256 and T — 1024 trading days 
correspond to values of the parameter s = mb = 0.04, 0.16, 0.64 and s = 2.56. 
We obtain that 

s T = m T b T w T/(400days) . (25) 

This determination of the time scale unfortunately is not very accurate and 
depends among other things on our assumption that all stocks are characterized 
by the same Tot (see below). The approach to a Gaussian fixed point of the 
negative tails is even slower 1 . The fact that we are approximating skewed 
distributions by symmetric ones may partly explain some of the discrepancies 
visible in Fig. 2. 

Systematic deviations also arise if the individual stock distributions scale 
less than perfectly for large T. If the parameter tot depends on the stock, 
the scaling among companies is broken on large time scales and the " lightest" 
stocks (the ones with the smallest tot) in this case dominate the average 
of Eq. (3) at large normalized returns. The observed flattening out of the tails 
of the averaged cumulative distributions for T > 16 days in this case could be 
due to some "lighter than average" stocks. 
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The overall quality of this one-parameter fit down to observed cumula- 
tive probabilities of 10 -5 is encouraging. The fact that st is approximately 
proportional to the trading time T furthermore suggests that the effective 
Hamiltonian could be time independent. For sufficiently long time intervals 
one would hope this to be the case for any parameterization of the evolution - 
most monotonically increasing quantities that might be relevant for the evolu- 
tion of a stock's price eventually do become proportional to the physical- (and 
to the trading-) time. 



5 Quantum Interference and Short Term Risk and Return 

For long time horizons T one expects that the mean return [It and the variance 
v T ~ hx/ftm?) of the returns are both proportional to T. We found that the 
time dependence of the distribution fits the empirical one if s = ra^T is 
approximately proportional to T. Taken together, this would imply that on 
time scales of T ~ years, 

\i T cx b T cx T, and m T = m <~ const. . (26) 

If one assumes time homogeneity of the returns, the quantum model in 
fact completely specifies the T-dependence of the parameters [it, br and tot 
for all T. Assuming that the distribution of historic returns of a stock over a 
time interval T does not depend on the initial time t is not very reasonable 
for single companies, but perhaps is tenable for short time intervals T once 
companies are grouped by their market capitalization - the procedure followed 
in rcf.[l]. The conditional probability density Pt{v\x) that the stock has a 
return y — x in the time interval T in this case is, 

Pr(y\x)\ t = p(y,t + T\x,t) cx \G{y,t + T;x,t)\ 2 , (27) 

where G(y,t+T; x, t) is the transition amplitude that a particular share traded 
at time t for a log-price x is traded at time t' = t+T > t for a log-price between 
y and y + dy. [The proportionality constant in general depends on T and t 
and is obtained from 1 = J p(y, t + T\x, t) dy.} 

The transition amplitude in the present model is* , 



G(y,t + T;x,t) =(y\Te- l l 6 ^^\x) 



— exp i(x - y)q + inr(t)q - b T (t)^q 2 + m T (t) 



.(28) 



tThe states |x) and \y) are normalized eigenvectors of the hermitian operator x that cor- 
responds to the log-price. The scalar product 17 with the previously introduced state \q) is 
(q\x) = e^ x . 
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An effective Hamiltonian that is diagonal in the wave number implies that 
G{y,t + T;x,t) and pr(y\x) depend on the return y — x only. If the returns 
were homogeneous in time, G(y, t + T;x, t) would not depend on the initial 
time t. This is an oversimplification of the dynamics and we will consider the 
possibility that the parameters in Eq. (28) vary slowly with the initial time, 
i.e. that the returns are homogeneous on time scales of the horizon T. 

The definition of Eq. (28) and the completeness of the basis of |j/)-states 
gives Trotter's formula for the transition amplitudes, 

G(z,T 1+ T 2 + t;x,t) = WTe-'ft™ 6 *^-'^*"^) 

J —oo 

/oo 
dy G(z, T 1 +T 2 + t; y, T\ + t)G(y, Ti + t; x, t) . (29) 

Eq. (29) is the mathematical expression of Huygens' principle describing the 
propagation of waves. Using Eq. (28) and performing the integral over the 
intermediate log-price y, one finds that for Eq. (29) holds only if, 

MT 1+ T 2 (t) = MT, (t) + ^T 2 {t + Ti) 

b Tl +T 2 (t) = bT 1 (t) + bT 2 {t + T 1 ) 
m Tl+T2 (t) - m Tl (t) = m T2 (t + Ti) = m(t) . (30) 

If the parameters depend little on the initial time (at least on time scales T± of 
interest), ^t 2 - [it 2 (t) and br 2 (*+Ti) - 6t 2 {t). The solution to Eq. (30) 

in this case is that ^r(^) and br(t) are both proportional to the time interval 
T, 

Hr(t) - Tjx{t), b T (t) ~ Tb(t), m T (t) = m(t) . (31) 

For notational clarity, the possible dependence of parameters on the initial 
time t will again be suppressed in the following. 

On a time scale of years, the T-dependence of the parameters is consistent 
with Eq. (26) and implies that both, the mean and the variance v\ ~ Tb/(2m) 
of long-term returns is proportional to T. On long time scales, Sharpe's ratio 
(r is the risk- free return rate), 

A T = MT - rT oc T 1 ' 2 for Tmb > 1 (32) 

VT 

is a monotonically increasing function of the investment horizon that is similar 
to that predicted by an iid process. On short time scales, the Sharpe ratio 
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of the present model does not increase as fast as for an iid process. Since 
the standard deviation vt ~ Tb is itself proportional to T, Sharpe's ratio in 
particular does not vanish for T —* 0, 

A =lim^^ = ^>0. (33) 

T^O V T b 

If returns were described by an iid process, short-term trading strategies would 
be extremely risky and should not find rational investors. The persistence of 
day-trading and of " excessive" trading by institutional investors 18 suggests that 
short-term strategies may not be prohibitively risky. The prolonged existence 
of such phenomena is more plausible if the specific risk, 1/At, remains finite 
for T — > 0. Although the risk per unit of excess return does rise, short- 
term investments are not penalized excessively by the price dynamics of the 
quantum model. [Short term investments nevertheless can be very expensive 
due to transaction costs.] 



Sharpe 's Ratio 




200 400 600 800 1000 1200 1400 
period in days 

Fig. 4: The time dependence of Sharpe's ratio. The points are from ref.[19] and give 
the observed Sharpe ratio of the Dow- Jones Industrials index. The solid line is the 
Sharpe ratio the quantum model predicts for an individual stock. The time scale 
has been calibrated by Eq. (25) and corresponds to the one of Fig. 2. The overall 
normalization of the Sharpe ratio of the stock was adjusted to the same order of 
magnitude as the index. For comparison, the dashed line gives the time dependence 
A(T) oc \pT of the Sharpe ratio that an iid process would predict. 

The Sharpe ratio of the historic returns on the Dow- Jones Industrials index 
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(DJI) in fact is fairly constant for periods up to T < 100 days 19 (and perhaps 
even declines somewhat). The results from Peters' analysis of the DJI index are 
shown in Fig. 4 and compared with the time dependence of the Sharpe ratio of 
a single representative stock in the quantum model. Although Peters studied 
the Sharpe ratio of an index rather than of individual stocks, his findings 
do indicate that the specific risk does not become prohibitively large as the 
investment period is shortened. 

To better see what a constant Sharpe ratio may mean to investors, compare 
the following two short-term investment strategies (0 < a < 1): 

51 An amount aW is invested in asset A and the remainder, (1 — a)W, is 
invested in asset B over a period T. 

52 The full amount W is invested in asset A for a time aT and for the 
remaining time, (1 — a)T, the full amount is reinvested in asset B. 

If the mean and the standard deviation of the returns on the individual assets 
A and B are uncorrelated and proportional to the investment period, both 
strategies carry the same short-term risk (as measured by the standard devi- 
ation of the overall returns) and have the same expected overall return. The 
equivalence of the two strategies does not depend on the ratio a nor on the 
standard deviations and mean returns of the individual assets**. If the Sharpe 
ratio approaches a finite value for T — > 0, a short-term investor therefore could 
be (almost) indifferent to choosing asset diversification (strategy SI) or time 
diversification (strategy S2). Since time diversification offers the option of 
choosing the reinvestment time at a later point in time, it will often be pre- 
ferred for short-term investments. [Note that the transaction costs incurred 
by both strategies could be similar and that time diversification may only be 
penalized by the fact that A > 0.] 

Less direct evidence for br oc T is provided by the probability density 
at zero return, pr(0) — pr{x\x). Eq. (15) implies that for T < 2 weeks, 
Pt(0) ~ in this model. With vt ~ &t oc T, the probability density at zero 
return on short time scales therefore is inversely proportional to the period, 
Pt(0) oc T -1 . vt oc \/T on the other hand would imply pr(0) oc T -1 / 2 . 
Neither of these power laws corresponds very well to pf, &P500 (0) oc T~ o n±om 
observed for the probability density at zero return of the S&P500 index 20 . But 
the stocks in an index are not perfectly correlated. Diversification decreases 
the variance of the index relative to the variance of the single-name stocks. 
Since correlations tend to decrease over time, the probability density at zero 

"The overall return and risk of either strategy docs of course depend on the returns on the 
individual assets and in general is optimal for a particular value of a only. 
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return on the index therefore increases relative to that of a single-name stock. 
A slower decay over time of the probability density at zero return for the index 
than for a single stock thus is more plausible than a faster one^. The prediction 
of the present model that Pt(0) ex T _1 for an individual stock thus may well be 
compatible with a probability density at zero return p!p P500 (0) cx T~ a - 71±0 - 03 
observed for the S&P500. 

One should emphasize that with the definition Eq. (27), Trotter's rela- 
tion Eq. (29) in general is not compatible with the relation for the transition 
probability densities of a Markov chain, 



because the absolute square of a sum of terms in general is not the sum of 
absolute squares of the individual terms. The difference is known as quantum 
interference. A Gaussian transition amplitude (and consequently a Gaussian 
transition probability density) is an important exception where Eq. (29) and 
Eq. (34) both hold. The absence of interference effects characterizes the (in- 
coherent) classical limit of a quantum system. Applying this criterion to the 
observed return distributions, the classical description of equity markets only 
becomes accurate on a timescale of years. 

Note that Eq. (34) may be violated in quantum theory even if the (quan- 
tum) process is mcmoryless and the Markov property holds for the condi- 
tional probability densities, i.e. if p(z,t"\y,t';x,t; . . .) = p(z,t"\y,t') for all 
t" > t' > t. This does not lead to any logical contradiction if price measure- 
ments at intermediate times have a non-negligible effect on the distribution. 
In a financial context this would means that trading a share to determine its 
price at an intermediate time results in a different distribution for the final 
outcome than if it were not traded. 

Although such back-coupling effects are well known and observable on 
financial markets, it is generally difficult to quantitatively include them in 
stochastic models. The quantum framework incorporates " measurement" - 
effects in a consistent and tractable (but perhaps rather narrow) fashion. There 
may be no need to model the impact from trading separately. 

6 Summary 

The observed 1 cumulative distributions of the short-term returns on single- 
name stocks are accurately modelled by pdf's of the form in Eq. (15). These 

ttThe probability density at zero return decays equally fast (oc l/y/T) for the stock and the 
index in the (uncorrelated) Gaussian case 




(34) 
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densities are the square of an amplitude that is an analytic function of the 
returns apart from two complex conjugate poles of opposite strength. The 
poles are located at z — /it + ivr and z* = jiT — »t, where \it is the mean 
return and vt is the standard deviation of the returns over the time period T < 
2 weeks**. This mathematically rather concise description trivially accounts 
for the observed 1 power law tails (with an exponent a = 3) and also implies 
scaling between the distributions of different individual companies. 

The fact that the distribution of short-term returns with finite variance 
approximately retains its shape over time scales T < 2 weeks rules out an iid 
process and is thus challenging to describe stochastically. The temporal evolu- 
tion of the corresponding amplitude on the other hand was found to be given by 
the remarkably simple linear operator H c s(t) of Eq. (18) that is proportional 
to the wave-number q. [The temporal evolution of the corresponding pdf (the 
absolute square of the amplitude) is not linear.] Since a dependence of the 
effective Hamiltonian on fractional powers of the wave number is unusual for 
quantum systems, a = 3 in this sense is a "natural" exponent for distributions 
with stable power law tails. 

The proposed quantum model gives a surprisingly quantitative and trans- 
parent explanation of the observed shape and temporal stability of the ob- 
served stock price fluctuations on short time scales. As Fig. 1 shows, the 
model quantitatively reproduces the average cumulative distribution observed 
for normalized 5 minute and daily returns over 7 orders of magnitude with no 
free parameters. 

By regulating the cusp of \q\ at q = and thus modifying the non-hermitian 
part of the Hamiltonian to a quadratic function of q near q = 0, the shape sim- 
ilarity of the evolution is broken and eventual convergence to a Gaussian prob- 
ability density over time horizons of several months to years is assured. With 
the modification of Eq. (19) the model qualitatively reproduces the observed 
distributions for the normalized returns 1 down to cumulative probabilities of 
1CT 4 for time periods up to 4 trading years. As discussed in section 4, the 
deterioration in the quality of the fit on long time scales may partly be due 
to the skewness of the observed distributions, which is more severe for longer 
time periods. It could also be due to some "lighter-than-average" stocks if 
the scaling of the individual stock distributions for T > 16 days is less than 
perfect. Such stocks would dominate the average of Eq. (3) at high normalized 
returns and explain the flattening out of the observed tails. 

Comparing the prediction of this model for the " decay-" , or trading- rate 
of a share with the observed 15 dependence of the trading rate on the imbalance 

^It is tempting to associate these singularities of the amplitude with long and short positions 
in the stock 
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between buyer- and seller-initiated trades, suggest that the imbalance number 
of ref.[15] is proportional to the wave- number. The imbalance in this case 
is the Fourier conjugate variable to the return xt (up to an overall scale). 
The quantum model then predicts a definite form for the distribution of the 
imbalance number - it in particular should have exponential, rather than power 
law, tails. 

By including sufficiently many "hidden variables" , that is factors, stochas- 
tic models may be able to reproduce the historic observations with similar 
accuracy. Indeed, the (possible) existence of hidden variables was at the heart 
of the early critique of quantum theory. For physical quantum phenomena, 
this alternative explanation has only recently been rejected by experimental 
verification of Bell's inequalities 21 . 

It is quite impossible to perform similar high-precision experiments on 
equities in order to reject the existence of "hidden variables". There in fact is 
no need to because it is almost self-evident that the return on a stock depends 
on many factors that have not been modelled. The question nevertheless is 
not just one of having a more efficient description of the dynamics. Unlike 
hidden variables describing physical phenomena, the factors that influence the 
dynamics of a stock are expected to change over time. It furthermore is not 
clear how economic factors like the gross national product influence the value of 
any given stock at any given point in time. The observed scaling of the return 
distributions for various stocks in different economic environments strongly 
suggests that all these " hidden" factors find their expression in the mean return 
and the variance of the returns. The dynamics that determines the shape of 
the return distributions on the other hand must be self-consistent and largely 
immune to the influence of " hidden" variables that are specific to a company 
and the economic and political climate. 

The conceptual advantage of a quantum description is that this consis- 
tent framework describes the observed return fluctuations rather well without 
separately modelling all the possibly influential, time-dependent, and hard-to- 
measure factors. Contrary to a purely descriptive "hidden variable" model, 
the consistency of a quantum model leads to a number of predictions. Their 
verification or falsification will ultimately decide whether this is a useful ap- 
proach. 
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A A Subordinated Stochastic Process for the 5 minute Returns 



In section 2 the T = 5 minute returns of a stock were found to be well described 
by the pdf of Eq. (15) with a — and standard deviation b — UT=5min (without 
loss in generality, we consider only the case of vanishing mean return). To 
characterize the directing process of a subordinated stochastic process 22 that 
would give this pdf, consider the following integral representation, 

%3 ^TxdXe-^ (35) 



7T(X 2 + 6 2 ) 2 7T 



00 e 1/(2t) (ft ex P 
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(36) 



One thus can interpret the pdf for the returns as the result of a random walk 
with variance tv"^, where t is itself a random variable drawn from the positive 
distribution 

e V(2t) 

" 7S» ' (37) 

Clark interprets the stochastic process r(t) whose increments are drawn from 
the distribution q(t) as "operational" time 9 . Note that q(t) has unit mean and 
infinite variance and that the random variable 1 jt is drawn from a T distribu- 
tion. In fact, the probability for a large variance, respectively operational time 
interval, itself falls off as a power law with an exponent of —(3/2). 

Although the subordinate stochastic process described here reproduces 
the observed distribution of returns for T = 5min by construction, it does 
not describe the observed time dependence of this distribution: the temporal 
evolution is not shape stable and a Gaussian fixed point is approached rather 
quickly. 
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